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We present a systematic study of the morphology of homoepitaxial InP films grown by metalor- 
ganic vapor-phase epitaxy which are imaged with ex situ atomic force microscopy. These films 
show a dramatic range of different surface morphologies as a function of the growth conditions and 
substrate (growth temperature, V/III ratio, and miscut angle < 0.6° and orientation toward A or 
B sites), ranging from stable step flow to previously unreported strong step bunching, over 10 nm 
in height. These observations suggest a window of growth parameters for optimal quality epitaxial 
layers. We also present a theoretical model for these growth modes that takes account of deposi- 
tion, diffusion, and dissociation of molecular precursors, and the diffusion and step incorporation of 
atoms released by the precursors. The experimental conditions for step flow and step bunching are 
reproduced by this model, with the step bunching instability caused by the difference in molecular 
dissociation from above and below step edges, as was discussed previously for GaAs (001). 

PACS numbers: 68.55.-a, 68.65.-k, 81.05.Ea, 81.10.Aj 



I. INTRODUCTION 

Homoepitaxial films of InP are commonly grown as 
buffer and cladding/waveguiding layers, setting the base 
for numerous types of semiconductor devices. 1,2 The mor- 
phology of these films can have a crucial impact on the 
quality of the overgrown material. For example, poor 
quality buffers could result in structural defects propa- 
gating to the device layers, reducing the carrier mobility 
or affecting optical properties, and thereby degrading the 
performance of the device. 3 ' 4 

In general, an epitaxial process is required to fabri- 
cate materials with reproducible properties of the highest 
quality in terms of crystal structure, purity, alloy com- 
position, surface morphology, etc. In most cases, mate- 
rial quality is controlled by the adjustment of the growth 
conditions, such as substrate temperature, gas/molecular 
flow and their ratios which, in metalorganic vapor-phase 
epitaxy (MOVPE), will affect precursors diffusion, sur- 
face decomposition processes, nucleation, diffusion of 
adatoms and their insertion at specific crystallographic 
sites and steps. For molecular-beam epitaxy, a well- 
founded conceptual and computational framework has 
emerged 6 since the pioneering work of Burton, Cabr- 
era, and Frank (BCF). 7 Systematic studies based on 
surface diffraction measurements and scanning tunnel- 
ing microscopy, in conjunction with extensive computa- 
tional modeling have identified and characterized many 
of the atomistic processes, even in the complex setting 
of III-V systems. 8-13 The corresponding development for 
MOVPE has been comparatively slow, due in part to 
the limited availability of in situ measurements. 14 Never- 
theless, systematic ex situ measurements of growth mor- 
phologies, again in conjunction with theoretical model- 
ing, have revealed the importance of precursor diffusion 



and decomposition at step edges in determining the na- 
ture and scale of surface morphologies. 15 

Epitaxy-ready substrates are prepared by cutting in- 
dividual slices out of the bulk crystal at the desired an- 
gle with respect to the main crystallographic plane. 16 In 
the case of InP, the crystallographic steps exposed after 
this cutting may be ideally terminated either with in- 
dium (so called A-stcps, with the substrate misorientcd 
toward the [1H]A planes), or phosphorus atoms (B-steps, 
with the substrate misoriented toward the [111]B planes), 
or any combination of these two. 17-19 Depending on the 
miscut, the exposed surface can be close to singular or 
contain a dense periodic array of crystallographic steps. 
Epitaxial overgrowth on such substrates may then pro- 
ceed (when truly "epitaxial" ) in one of three basic growth 
modes: step flow, where the overgrowing layer is advanc- 
ing each exposed step at the same rate, creating an exact 
copy of the underlying surface; step bunching, where ter- 
race formation is observed due to clustering of individual 
steps; or island formation, when the growth of new layer 
is initiated only at the step edges, but also between steps, 
in which case no long range order is observed. If the 
growth conditions are not optimal, due to defect forma- 
tion, the surface morphology can be corrupted and will 
not resemble any of the three main modes. For most ap- 
plications, step flow is the desired mode, as it leads to cre- 
ation of atomically flat layers, without irregularities (such 
as, in an extreme example, antiphase boundaries). 20 

On the other hand, the role of small substrate mis- 
orientations for the III-V system have been recently re- 
vealed for arsenide alloys (e.g. GaAs, AlInAs or InAs) 
grown by MOVPE, 15 ' 21 ~ 24 where complex surface orga- 
nization, striking effects and state-of-the-art results on 
the material optical/transport properties have been re- 
ported. In those reports, small changes of growth con- 
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ditions significantly affected the epitaxial morphology. 
As usual, the growth conditions revolved around a few 
parameters, where the common variables were substrate 
miscut, growth temperature, growth rate and molar ra- 
tio between precursors injected into the reactor chamber 
(reactor pressure, carrier gas and precursor choice were 
fixed in the given setup). Crucially, the authors found a 
strong correlation between surface organization and ma- 
terial properties. What is even more relevant is that the 
authors linked the unusual variety of surface organization 
(e.g. from islanding, through step flow to periodic step 
bunches just by varying substrate miscut) , to the intrinsic 
two-step MOVPE growth process: first adsorption, dif- 
fusion and decomposition of molecular precursors at step 
edges and, subsequently, adatom diffusion and incorpo- 
ration. Their model reproduced the condition for island 
nucleation as well as for the step bunching instabilities, 
linking this last process to the difference in molecular 
dissociation from above and below step edges. 15 

Despite the broad and primary commercial interest 
and exploitation of InP, 25 only limited data are avail- 
able in the literature on InP grown by MOVPE. Only 
a small range of parameters/surface organization appear 
to have been correlated. For instance, the influence of 
small differences of substrate misorientation was partly 
discussed in Ref. [26] (for substrates with offcuts of 0.02°- 
0.25° toward (Oil)) and Ref. [27] (for singular wafers 
and with misoricntations of 0.2° toward [110], [111] A 
and [111]B). Morphologies created during homoepitaxy 
on singular wafers (with misorientation below 0.1°) were 
also investigated in Ref. 28 with respect to the direction 
of the miscut. In Ref. [28], the authors analyzed the ef- 
fect of temperature on step bunching on vicinal surfaces 
in a larger miscut range (from 0.2° to 2° toward [111]A). 
Nevertheless, there is no comprehensive study of the or- 
ganization of the MOVPE-induccd crystallographic steps 
on vicinal surfaces, so a fundamental understanding of 
the fundamental kinetic processes is still lacking. 

In this paper, we report a systematic study of the im- 
pact of growth conditions on the surface morphology of 
InP surfaces with small misoricntations during MOVPE. 
We show several surprising effects and a remarkably 
broad variety of surface organizational modes, well be- 
yond what has been observed in the arsenide (GaAs) 
counterpart. The organization of our paper is as follows. 
The procedure used for our experiments is described in 
Sec. II. The surface morphologies observed on vicinal 
InP(OOl) are described in Sec. Ill, with the presentation 
of the results grouped according to growth rate and V/III 
ratio (Sees. Ill A, IIIB, IIIC) and the effect of doping 
(Sec. HID). These results are analyzed in Sec. IV, which 
includes a theoretical discussion based on an extension 
of the model in Ref. [15]. One of the interesting aspects 
of our experiments is that we can examine the effect of 
the group-V species on the various steps of the growth 
kinetics. We summarize our results in Sec. V, where we 
also discuss future modelling strategies. 



II. EXPERIMENTAL PROCEDURE 

All samples in this study were grown by MOVPE at 
low pressure (80 mbar) in a commercial horizontal reactor 
with purified N2 as the carrier gas and trimethylindium 
(TMIn) and phosphine (PH3) as precursors. 29 300-nm- 
thick single layers of InP were grown on perfectly ori- 
ented or slightly misoriented InP (001) substrates, which 
in most cases were semi-insulating Fc-dopcd, with ref- 
erence samples grown on Zn- (p-typc) and Si-doped (n- 
type) wafers (we anticipate no observed differences to 
result from the substrate doping, as no diffusion to the 
epitaxial layer is expected). The majority of the InP lay- 
ers were grown as nominally undoped. The unintentional 
doping level is estimated to be well below 5 x 10 15 cm~ 3 , 
which is our detection threshold. We have also grown 
several samples using dicthylzinc (DEZn) and disilanc 
(Si2Hg) as precursors to investigate the impact they have 
on the surface morphology when incorporated into the 
InP matrix. Doping levels were 5 x 10 17 cm~ 3 for Zn and 
9 x 10 17 cm~ 3 for Si and constant across the layer, as 
confirmed by electrochemical capacitance voltage mea- 
surements. Such carrier concentrations are considered 
moderate for InP-bascd semiconductor devices. 30 

Growth conditions were varied systematically across 
samples: the estimated actual growth temperature Tq 
varied in the range 520°C to 720°C, the V/III ratio R 
from 150 to 450, and sample misorientation from nom- 
inally "perfectly oriented" to 0.6° (all with a tolerance 
of ±0.02°) toward [111] A or B. The growth rate G was 
kept constant at 0.7 /zm/hr for most of the samples, with 
reference/comparison growth carried out at 0.35 /Ltm/hr 
and 1.4 /im/hr. All together, more than 150 samples were 
grown in various conditions for this study, with several 
control growth runs performed to ensure the full repro- 
ducibility of the reported results. In all cases, partic- 
ular attention was paid to reactor environment quality 
and to temperature control by growing the samples only 
with reactor walls already baked and covered by previous 
growth runs. The growth temperature was estimated by 
cmissivity corrected pyrometry. The relevant details are 
referenced below when a given example is discussed. 

All epitaxial growth resulted in mirror-like surfaces. 
The samples were first inspected with an optical mi- 
croscope in (Nomarski) differential interference contrast 
(NDIC), as well as in dark-field mode. The vertical reso- 
lution of NDIC is limited and a subsequent detailed mor- 
phological study was performed with atomic force mi- 
croscopy (AFM) in tapping/non contact mode at room 
temperature and in air. The results are presented in vari- 
ous scan sizes to emphasize the features relevant to a par- 
ticular growth morphology. While we always measured 
our samples at a short interval from removal from the 
reactor, the surface features did not appear do degrade 
significantly over several weeks, an indication that, as for 
GaAs, the surface oxide develops slowly enough and con- 
formally to the original structure in the first months of air 
exposure. The quality/crystallinity of the grown material 
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TABLE I. The growth rate, V/III ratios, misorientation, and growth temperature of InP substrates produced during homoepi- 
taxy by MOVPE. The samples are labeled for each set of growth conditions, together with an abbreviation for the observed 
morphology: discrete islands (I), diffuse islands (DI), step flow (SF), "normal" step bunching (SB), cliffs (C), braids (B), and 
defects (D). 



Growth rate 


V/III ratio 


Substrate 






Growth temperat 


ure 






520°C 


560° C 


585° C 


610°C 


630 C 


655°C 


685° C 


720° C 






"perfectly oriented" 


1 A 

1A 

I 


2 A 

I 


3A 
I 


4A 

SF 


5A 

SB 


a A 
() A 

SB 


T A 

7A 

SB 


O A 
SA 

D 






0.2° toward [111]A 


IB 

DI 


>, i > 

ZD 

DI 


3B 

C 


4B 

C 


5B 

C 


1. 1 1 

B 


7B 

D 


,i 1 1 
aB 

D 




150 


0.4° toward [111]A 


1C 

DI 


2C 
DI 


3C 

C 


4C 

C 


5C 

C 


B 


D 


D 






0.4° toward [111]B 


ID 

DI 


.lit 
2L> 

DI 


3D 

C 


4D 

C 


5D 

B 


oD 

D 


— 1 1 

7D 

D 


SD 

D 


0.7 fxra/hr 




0.6° toward [111]A 


IE 

DI 


2E 

DI 


3E 

C 


4E 

C 


5E 

B 


6E 

SB 


7E 

D 


8E 
D 




"perfectly oriented" 


HA 

SF 




10A 

SF 




11A 

SB 


1.1 A 

12 A 

SB 




1 O A 

SF 






0.2° toward [111]A 


9B 

SB 




10B 

B 




11B 

B 


12B 

B 




13B 

SB 




450 


0.4° toward [111]A 


9C 

SB 




IOC 

B 




11C 

B 


12C 

B 




13C 

SB 






0.4° toward [111]B 


9D 

SF 




10D 

SF 




11D 

SB 


12D 

SB 




13D 

SF 






0.6° toward [111]A 


9E 

SB 




10E 

B 




HE 

B 


12E 

B 




13E 

SB 



TABLE II. The growth rate, V/III ratios and misorientation 
of InP substrates for homoepitaxial growth at Tq = 630° C 
during MOVPE. The samples are labeled for each set of 
growth conditions, together with an abbreviation for the ob- 
served morphology used in Table I. 



Substrate 


G = 0.35 fim/hi 


G = 1.4 /xm/hr 


i? = 450 


R = 150 


i? = 450 


"perfectly oriented" 


14A 

SF 


15A 

SB 


16A 

SF 


0.2° toward [111]A 


14B 

B 


15B 

B 


16B 

SB 


0.4° toward [111]A 


14C 

B 


15C 

B 


16C 

SB 


0.4° toward [111]B 


14D 

SF 


15D 

B 


16D 

SF 


0.6° toward [111]A 


14E 

SB 


15E 

B 


16E 

SB 



was also confirmed by X-ray diffraction measurements. 



III. RESULTS 

Our experiments revealed a surprising variety of InP 
surface morphological organization during homocpitaxy 
by MOVPE. The growth parameters, which are the 
growth rate G, V/III ratio R, misorientation angle and 



TABLE III. Summary of surface morphologies of the samples 
in Tables I and II. 



Surface morphology 


Sample number 


Islands 


Discrete (I) 


1A, 2A, 3A 


Diffuse (DI) 


IB, 1C, ID, IE, 2B, 
2C, 2D, 2E 


"Plain" 


Step flow (SF) 


4A, 9A, 9D, 10A, 10D, 
13A, 13D, 14A, 14D, 16A, 
16D 


"Normal" step 
bunching (SB) 


5A, 6A, 6E, 7A, 9B, 
9C, 9E, 11A, 11D, 12A, 
12D, 13B, 13C, 13E, 14E, 
15A, 16B, 16C, 16C, 16E 


Step bunching 
with 
long-range 
periodicity 


Cliffs (C) 


3B, 3C, 3D, 3E, 4B, 
4C, 4D, 4E, 5B, 5C 


Braids (B) 


5D, 5E, 6B, 6C, 10B, 
10C, 10E, 11B, 11C, HE, 
12B, 12C, 12E, 14B, 14C, 
15B, 15C, 15D, 15E 


Defects 


Defects (D) 


6D, 7B, 7C, 7D, 7E, 
8A, 8B, 8C, 8D, 8E 



direction, and growth temperature Tq for each sample 
are compiled in Tables I and II, with labels that will be 
used in the following discussion. The samples associated 
with each type of morphology are listed in Table III. We 
first discuss the results for the growth of nominally un- 
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FIG. 1. (Color online) AFM images (signal amplitudes) of the 
top surface of samples grown with R — 150, at Tq = 520° C 
(left panel) and Tq = 560° C (right panel), on substrates that 
are perfectly oriented (±0.02°) (1A, 2A), misoriented by 0.2° 
toward [111] A (IB, 2B) and by 0.4° toward [111]B (ID, 2D). 



doped layers for different growth rates and V/III ratios, 
before considering the effect of doping. 



A. G = 0.7 /xm/hr, R = 150 

On samples 1A-3A (low Tq, no intentional misorienta- 
tion) , we observed the random nucleation of distinct and 
separated epitaxial islands between the crystallographic 
steps (when discernible) of the substrate (Fig. 1). This 
morphology will be referred to as "islanding" and signi- 
fied by I. On samples IB-IE and 2B-2E (low Tq, with 
intentional substrate misorientation) , however, a disor- 
ganized mixture of step flow and step bunching devel- 
oped on which the steps of the original substrate could 
not be easily discerned (Fig. 1). This morphology, called 
"diffuse islanding", will be signified by DI. This type of 
morphology is similar to that observed in models and 
scans with scanning tunnelling of growth on misoriented 




FIG. 2. (Color online) AFM images (signal amplitudes in top 
row and corresponding three-dimensional (3D) height recon- 
structions in bottom row) of samples showing cliff-like step- 
bunching. Growth was performed on substrates with a mis- 
orientation of 0.4° toward [111]A (4c), 0.4° toward [111]B 
(4D), and 0.2° toward [111]A (5B) at growth conditions with 
R = 150, G = 0.7 jum/hr, T G = 610°C (4C, 4D) and 
T G = 630° C (5B). 



surfaces during MBE 8,31 in a transition regime between 
island nucleation and growth and step flow. When the 
growth temperature was raised we observed more orga- 
nized surface features. The samples grown on perfectly 
oriented wafers showed morphologies close to step flow, 
but with some residual islands (4A), step bunching of 
several neighboring steps without significant long range 
organization (5A), or in close grouping of two subsequent 
steps (6A.) The step-bunched morphology will be de- 
noted as SB. 

On samples with an intentional misorientation, signif- 
icant step bunching was observed at lower temperatures 
than for the perfectly oriented surfaces. Samples 3B-3E, 
4B-4E and 5B-5C showed the formation of "cliffs", de- 
noted as C, made of up to steep 40 monolayer stacks 
(corresponding to 10- nm- high edges). The general trend 
is similar for all these samples, but the observed features 
were the most distinct on samples 3D, 4B-4E and 5B 
(Fig. 2). When investigated in detail, the cliff morphol- 
ogy shows a very dense succession of nm step bunches, 
rather than atomically sharp edges. 

Different morphologies were obtained for Tq > 630° C. 
This was first evident for samples with largest misori- 
entation (5D, 5E), and then at Tq = 655°C for all A- 
misoriented samples (6B, 6C, 6E), while B-misoricntcd 
(6D) and perfectly-oriented (6A) samples showed moder- 
ate step bunching, but with an indication of defect forma- 
tion. Strongly step-bunched surfaces were still observed, 
but the bunched steps were not simply overlapping; in- 
stead, while maintaining a closely grouped and periodic 
structure, they showed significant numbers of single steps 
between the bunches, creating a more "braid- like" picture 
on a flattened image (Fig. 3). Moreover, on a scale of 10- 
30 /im, the "braids" were seen to dissolve into single steps 
which subsequently plaited into the next "braid" ahead. 
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FIG. 3. (Color online) AFM images (signal amplitudes in top 
and bottom row and 3D height reconstruction in central row) 
of samples showing braid-like step bunching. The samples 
were grown on substrates with a misorientation of 0.4° (5D) 
and 0.6° toward [111] A (5E) with R = 150 and G = 0.7 ^m/hr 
at Tg = 630°C. The top and central rows contain images 
in close zoom, while the bottom row shows the large-scale 
organization of the features. 

This morphology will be signified by B. 

For the highest growth temperatures, we observed sig- 
nificant disruption of the sample surfaces. Indentations 
with a depth of up to 5 nm, corresponding to approx- 
imately 20 monolayers (MLs), were measured on sam- 
ples 8A-8E, with smaller dimples and step pinning ob- 
served also at slightly lower temperatures (sample 6A, 
6D and 7A-7E). Plan view images, lines scans, and three- 
dimensional images of such surfaces are shown in Fig. 4. 
These defected surfaces are signified by D. 

B. G = 0.7 /xm/hr, R = 450 

Samples grown under higher phosphine flow showed 
a much less diverse range of morphologies than those 
with lower flows. Starting from growth temperatures 
as low as 520°C we were able to identify organized step 
behavior on the surfaces of all samples. The singular 
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FIG. 4. (Color online) AFM images (signal amplitudes in 
top row, corresponding height profiles in central row and 3D 
height reconstructions in bottom row) of sample surfaces. Im- 
ages refer to samples grown at R = 150, G = 0.7 ^im/hr, 
T G = 655° C (6A, 6D) and T G = 720° C (8A), using substrates 
that are perfectly-oriented (6A, 8A) or misoriented by 0.4° 
toward [111]B (6D). Samples 6D and 8A show defected sam- 
ple surfaces, while the sample 6A shows step pinning (sharp 
change in the step edge line). 



and B-misoricnted samples (9A, 9D) showed step-flow, 
while the A-misoriented samples showed normal "disor- 
dered" step bunching (with step jumps of a few MLs, 
up to 2 nm) without long-range periodicity (9B-9C, 9E). 
Nevertheless, a slight increase in the growth temperature 
resulted in long-range step organization. While singu- 
lar and B-misoricnted samples (10A-12A and 10D-12D) 
seemed not to be affected, all of the A-misoriented sam- 
ples (10B-10C, 10E) showed significant step bunching, 
which became even more pronounced at higher temper- 
atures (11B-11C, HE, 12B-12C, 12E), resembling the 
morphologies in Figs. 2 and 3. We again observed the 
plaiting braid-like morphology, even with increased or- 
dering, which remained up to our highest growth tem- 
perature (T G = 720°C), as shown in Fig. 5. 

Contrary to what was previously observed (defected 
areas), in these growth conditions, the high growth tem- 
perature resulted in very smooth surfaces, with step flow 
or limited/normal step-bunching on high off-cut wafers 
(Fig. 6), which might be an indication that the previously 
observed indentation might be related to an insufficient 
supply of phosphorus to the surface. 

C. G = 1.4/0.35 /im/hr, R = 450/150 

The effect of varying the growth rate and the phos- 
phine flow rate was investigated only at Tq = 630° C, 
with the morphologies summarized in Table II. Notably, 
when the growth rate was increased, we observed an im- 
provement in the uniformity of the surface. The effects 
of the change (at G — 1.4 /im/hr compared to our stan- 




FIG. 5. (Color online) AFM images (signal amplitudes in top 
and bottom row and 3D height reconstruction in central row) 
of samples showing braid-like step bunching. These images 
refer to samples grown on substrates with an off-cut of 0.2° 
toward [111]A at R = 450, G = 0.7 fim/hx, T G = 630°C (11B) 
and Tq — 655°C (12B). The top and central rows contain 
images in close zoom, while the bottom row shows the large- 
scale organization of the features. 




FIG. 6. (Color online) AFM images (signal amplitudes) of 
the top surface of samples grown at 720° C and R — 450 on 
substrates that are perfectly oriented (13 A) and misoriented 
by 0.2° toward [111]A (13B) and 0.4° toward [111]B (13D). 



dard of G = 0.7 /um/hr) was consistent for samples with 
both high (R = 450) and low (R = 150) phosphine flow 
rates: the large step bunching was softened, resulting 
in the absence of long-range periodicity for the samples 
grown with R = 450, while a softening of the cliff-like 
edges on samples growth with R = 150 (Fig. 7). On 
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FIG. 7. (Color online) AFM images (signal amplitudes) of 
the top surface grown on samples misoriented by 0.4° toward 
[111]A at R = 450 (11C) and R = 150 (5C) at G = 0.7/xm/hr, 
and G = 1.4^m/hr (16C, 15C). Bunched areas in samples 5C 
and 15C are ~20 MLs in height with a shape change. 5C 
shows more sharp cliffs, while 15C is more slope-like. 



the other hand, reducing the growth rate to 0.35 ^tm/hr 
seemed to have no effect on the morphology compared to 
samples grown with 0.7 ^tm/hr. 

It is interesting to note that, as expected, the sam- 
ples showing island formation, step flow and "normal" 
step bunching appear to be featureless under the optical 
microscope. Nevertheless, the significantly step-bunched 
samples had surface morphologies which enabled char- 
acterization with either dark- field microscopy or N-DIC. 
While the N-DIC images clearly showed the lines cor- 
responding to the feature edges, the dark-field configu- 
ration accentuated the intertwining of the "braid-like" 
bunches (Fig. 8). This is significant for standard labora- 
tory practice, as often first quality control of the growth 
is performed under optical microscope. Thus it is impor- 
tant to note, that observed features do not necessarily 
correspond to defected layer and might be originating 
from surface morphology. 



D. Doping 

Depending on the application, the structural design 
of a device often calls for conductive InP layers. This 
is realized by the intentional doping of the bulk materi- 
als, often with elements such as zinc (Zn) or silicon (Si). 
In InP the Si could in principle introduce either a donor 
state by replacing In, or an acceptor state by replacing P, 
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FIG. 8. (Color online) Long-range organization of the surfaces 
of samples 4C and 6C imaged by AFM (signal amplitude, top 
row), N-DIC microscopy (middle row) and dark-field optical 
microscopy (bottom row). Optical and AFM images are of 
different scale, as indicated by the scale bars. 



but the donor state is more energetically favorable, and 
therefore dominant. 32 The incorporation of Si is propor- 
tional to the disilanc flow and the incorporation efficiency 
is high, i.e. the silane flow needed to achieve our target 
doping levels is an order of magnitude lower than the 
TMIn flow. 33 However, Zn doping normally has a very 
low incorporation efficiency, 34,35 so less than 0.2% of the 
molecules injected into the reactor were incorporated into 
the crystal matrix with a molar ratio TMIn/DcZn sw 50. 

In Table IV we list the growth parameters for the sam- 
ples discussed in this section. Figure 9 presents a compar- 
ison of the morphology of layers grown on semi-insulating 
substrates with those discussed in preceding sections that 
have the same misorientation. For this study we used 
R = 520 (slightly different from previously used values, 
as optimized for different experiments in our laboratory, 
nevertheless comparable to the high V/III ratios reported 
before), Tq = 630°C and G = 0.7 /xm/hr. Good surface 




FIG. 9. (Color online) AFM images (signal amplitudes) of the 
top surface of undoped (left panel), Si- and Zn-doped sam- 
ples (central and right panel, respectively) grown on perfectly 
oriented substrates (top row), and substrates misoriented by 
0.2° toward [111]A (middle row) and by 0.4° toward [111]B 
(bottom row). 



TABLE IV. Misorientation and layer dopant for InP sub- 
strates at To = 630° C and R = 520. The samples are labeled 
for each set of growth conditions, together with an abbrevia- 
tion for the observed morphology used in Table I. 



Substrate 


Dopant 


Zn 


Si 


"perfectly oriented" 


17A 

SF 


18A 

SB 


0.2° toward [111]A 


17B 

SF 


18B 

B 


0.4° toward [111]A 


17C 

SB 


18C 

B 


0.4° toward [111]B 


17D 

SB 


18D 

B 


0.6° toward [111]A 


17E 

SB 


18E 

B 



quality was obtained for all growth runs, as expected. 
However, differences were observed when step-step inter- 
actions are of interest. 

Remarkably, Zn-dopcd layers showed significantly less 
step-bunching, while the Si doping seemed to have less 
of an impact, even if it appeared to have intensified 
kink formation and affected growth on B-misoriented 
substrates more than that on A-misoriented substrates. 
While for undoped cpilayers with similar growth condi- 
tions (R = 450) we observed bunching of two steps on 



singular wafers, braid-like step bunching on A-type sub- 
strates, and normal step bunching on B-type substrates, 
the incorporation of Zn atoms led to ideal single step-flow 
on perfectly oriented substrates and normal/modest step 
bunching (of only a few monolayers) on all misoriented 
wafers, without any unusually large terrace formation. 
On the other hand, Si-doped layers showed "braid" -like 
step bunching, even on B-typc substrates, rarely seen in 
undoped layers (with the exception of sample 5D). The 
use of Zn- or Si-doped substrates (as we anticipated ear- 
lier in our contribution), seemed to have no effect on the 
surface morphology - the features were consistent with 
those seen on semi-insulating substrates (not shown). 



IV. DISCUSSION 
A. Terrace lengths and step bunches 

The wide range of growth conditions used in our in- 
vestigation has revealed several identifiable trends in the 
resulting surface morphologies. The V/III ratio had 
the strongest impact on the morphology: for growth at 
R = 150, the temperature window for obtaining good 
surface morphology was extremely narrow. Growth tem- 
peratures lower than 585°C were not sufficient to allow 
the formation of crystallographic steps while, at 655° C, 
we already observed the appearance of surface defects 
(Fig. 10). Increasing to R = 450 expanded the allowed 
range to essentially all temperatures we used - no defects 
were observed even at 720° C, nor were the surfaces pro- 
duced at 520°C significantly worse than those at higher 
temperatures (Fig. 10). We observed a lengthening of 
the terraces with temperature and, of course, a decrease 
of this length with larger substrate misorientation angle. 
(Figs. 11 and 12). An increase of the phosphine flow to 
R = 650 did not bring about any observable change in 
surface morphology. We also observed good results for 
growth with R = 370 (not shown). Cliff- like step bunch- 
ing was observed only for lower phosphine flows. 

The growth rate also has a significant impact on the 
surface features. The terrace formation observed on lay- 
ers grown at G = 0.35 /x/hr and 0.7 /im/hr was trans- 
formed into more conventional step bunching or step flow 
for G = 1.4/im/hr and R = 450, and the cliff-like kinks 
were "softened" into braids for R = 150 (Fig. 7). We 
nevertheless caution the reader that, as we have already 
mentioned, at high growth rates, we limited our studies 
to a specific growth temperature, and that further work 
is necessary to draw broader conclusions. 

The defects observed for growth at high temperatures 
could be easily associated with desorption of phosphine 
from the surface. What supports that conclusion is the 
fact that higher phosphine flow cured the sample surface, 
creating good morphologies even at very high growth 
temperatures. Additionally, the fact, that B-type sub- 
strates (containing phosphorous-terminated steps) have 
shown defect formation at lower temperatures than A- 
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FIG. 10. (Color online) Comparison of different morphologies 
as described in the main text. AFM images (signal ampli- 
tudes) of the top surface of samples grown with R = 150 (left 
panel) and R = 450 (right panel), at Tc = 585° C (top row), 
630° C (middle row) and 720° C (bottom row). 



type samples, seems to suggest the same picture. 

Some aspects of the differences observed in the growth 
on substrates with A- and B-steps can also be understood 
from the generic (i.e. not materials specific) properties of 
these steps. A-steps are group-III-terminated (In for InP) 
and B-steps are group- V-terminated (P for InP). These 
steps have very different formation energies, which is re- 
flected in their kinetics. The A-steps have a lower energy 
of formation than the B-steps. 36 and an In atom has a 
higher detachment barrier from a B-step than from an A- 
step. Hence, an In adatom can attach and detach from an 
A-step much more easily than from a B-step, so A-steps 
are smoother (straighter) than B-steps. 37 This seems to 
be supported by the images presented earlier, and espe- 
cially highlighted in Fig. 13 for surfaces misoriented by 
0.4° along the A and B directions. 
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FIG. 11. (Color online) Dependence of terrace lengths af- 
ter growth on the initial substrate misorientation angle and 
growth temperature for samples grown with R = 450 and 
G = 0.7 /Ltm/hr. 



B. Stable and unstable modes on misoriented 
surfaces 



We can obtain a qualitative understanding of the 
morphological variations on A-surfaces of misoriented 
InP(OOl) by appealing to an earlier study of growth on 
GaAs (001) by MOVPE. 15 The trend observed there on 
surfaces misoriented by 0° (±0.02°), 0.2°, and 0.4° was 
islanding, step flow, incipient step bunching, and well- 
developed step bunching, respectively. These results can 
be explained by a model that supposes that the decom- 
position of trimethylgallium is greater from above than 
below steps, which produces a net downhill current, as 
would an inverse Ehrlich-Schwocbel barrier. 38 ' 39 The ki- 
netics of the group-Ill adatoms were comparatively unim- 
portant under these growth conditions. However, the ex- 
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FIG. 12. (Color online) Comparison of measured distances 
between bunches on samples grown with R = 450 (circles) 
and R = 150 (filled squares) on substrates with the indicated 
misorientation angles. 



pcriments reported here have been carried out under a 
much wider range of growth conditions, including large 
variations in the phosphine flow, which will enable a 
much more extensive study of the growth kinetics. For 
this reason, we will retain the full range of kinetic pro- 
cesses for both precursors and adatoms, including asym- 
metric step-edge incorporation. 

Our model 15 is an extension of the basic BCF theory 7 
for the surface concentrations n{x, t) of group-Ill precur- 
sors (i.e. TMIn for InP) and c(x,t) of group-Ill adatoms 
(i.e. In) that includes the deposition rate of of precursors 
to the surface, their subsequent surface diffusion and pos- 
sible dcsorption from the substrate, their decomposition 
at step edges, the release of group-Ill atoms from their 
precursors and the subsequent adatom surface diffusion. 
The coupled reaction-diffusion equations for n and c that 
describe these processes are 



dn d 2 n n 

-TT7 = Dm^-o KTI + F , 

at ox z t 



dc p d 2 c 
dt dx 2 



(1) 
(2) 



where Dm and Da are, respectively, the surface diffusion 
constants of the precursor and the adatoms, K and r _1 
are the decomposition and dcsorption rates, respectively 
of the precursor on the terrace, and F is the effective de- 
position flux of the precursor. The quantity nn therefore 
acts as an effective spatially-dependent deposition flux 
for group-Ill atoms. Although these equations are ap- 
plicable for growth of a two-dimensional surface, in the 
interest of obtaining an analytic theory, we will confine 
our selves to a one-dimensional surface (Fig. 14). This, 
of course, neglects any step meandering which, as the 
images in the preceding section reveals, is an interesting 
issue in its own right. 

The reaction-diffusion equations (1) and (2) are sup- 
plemented by boundary conditions at the leading and 
trailing edges that bound each terrace. On the nth ter- 
race, for which x n < x < x n +i, the boundary conditions 
for the precursors are 

D M n x (x n ,t) = P M n{x n ,t) , (3) 

-D M n x (x n+ i,t) = /3 M n(x n+1 ,t) , (4) 

and for the adatoms, 

D A c x {x ni t) = /3^[c(x n ,t) - c ] , (5) 

-D A c x {x n+ll t) = pj[c(x n+1 ,t) - c ] , (6) 

where Co is the equilibrium adatom concentration at the 
step edge. The boundary conditions for the precursor 
stipulate that a molecule incident on step from above 
(+) or below (— ) decomposes and the group-Ill atomic 
constituent incorporated at a rate proportional to (3 M . 
The boundary conditions for the group-Ill adatoms state 
that atoms incident on a step from above or below are 
incorporated into the solid at a rate proportional to . 
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FIG. 13. (Color online) AFM images (signal amplitudes in the top row, 3D height reconstructions in bottom row) of samples 
grown with G = 0.7^m/hr at T G = 630°C and R = 150 (5C, 5D) or R = 450 (11C, 11D) on substrates misoriented by 0.4° 
toward [111]A (5C, 11C) or [111]B (5D, 11D). 




FIG. 14. The positions x n of steps on a one-dimensional sur- 
face and the surface currents from above ( — ) and below (+) 
each step that drive step motion and determine the stability 
of the step train. 



The morphological evolution of a misoriented surface is 
described by the changes in the positions x n of the steps 
as the result of absorbing material from surface currents 
(Fig. 14), 



J + {L n ) + J-(L n _i), 



(7) 



where L n = x n+ i — x n and = x n — x n -i are the 

lengths of the upper and lower terraces, respectively, at 
x n , and the surface currents J ± are obtained by solving 
(l)-(6) for each terrace in the system. The corresponding 
equations for the terrace lengths L n are 



dt 



~~dt df 

J+(L n+1 ) + J-(L n ) 



J + (L n )- J-(L„_i) 



(8) 



These equations have the stationary solution L n = L for 
all n, where L = aj tan# is the average terrace length 
determined by the misoricntation angle 9, where a is the 
step height which, for InP(OOl) is the height of an In-P 



bilayer. The central question in this paper is whether 
the regular train is stable to step bunching. This can 
studied most expediently with a linear stability analysis 
by considering small deviations A n from the regular step 
train: L n = L + X n . By retaining terms only to first order 
in the A„ in the Taylor series of J ± , 

J ± (L n±1 ) = J ± (L)+ dJ 



dL 



An 



±i 



(9) 



where the subscript on the derivatives indicates evalua- 
tion with the regular step train, the linearized form of 
(8) is obtained as 



dXn 

~df 



dJ+ 



dL 



— A n ) + 



dJ~ 



<JL 



(A n -A„_i). (10) 



With the step displacements expressed as Fourier modes, 

A» = itfce^-*"** , (11) 

the growth or decay rate Rk of the fcth mode, which is 
the real part of cuk , is 

i?, = 2sin 2 (iA : L)^ , (12) 
dL 

where J S {L) = J~{L) - J + {L). The stability of the reg- 
ular step train is determined by the sign of Rk which, in 
turn, is determined by the sign of J' S (L): if Rk < 0, the 
regular step train is stable, while if Rk > 0, it is unstable 
to step bunching. The most unstable mode corresponds 
to k = n/L, with eigenvector (1, — 1, 1, —1, . . .), so step 
bunching is initiated by the pairing of adjacent steps. 



C. Stability of step flow on misoriented InP(OOl) 

The decay rate Rk is determined from the station- 
ary solution of (l)-(6). The general solutions to (1) 
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and (2) are straightforward to determine, but the sub- 
sequent calculation of the currents and their derivatives 
are quite lengthy. Therefore, the generation and ma- 
nipulation of these solutions have been carried out in 
Mathematica. 40 The rates T n of all kinetic processes in 
our model have an Arrhenius form, T n = v n e~ En l kBT ', so 
determining Rk necessitates assigning values to the pref- 
actor v n and barrier E n for each process. As our model 
does not explicitly include the group- V kinetics, these pa- 
rameters have an implicit dependence on the phosphine 
flow rate. 

Our model will be used to analyze the morpholo- 
gies in Table I, as this provides the most comprehen- 
sive temperature- misorientation (toward [111] A) matrix 
of data. These data and our fits are shown in Fig. 15, 
with the Arrhenius parameters for the kinetic coefficients 
appearing in (l)-(6) compiled in Table V. The focus of 
most of our discussion will be on processes occurring at 
step edges, but the prefactor for the surface diffusion of 
TMI merits some comment. The value v = 10 14 , which 
may seem quite high in comparison with the range 10 12 - 
10 13 typically used in such theories and simulations, actu- 
ally conforms with estimates 41 for large molecules. Such 
large prefactors are due to jumps over several lattice cites, 
which is in contrast to the usual assumption of jumps of 
a single lattice unit. Otherwise, the precise values in Ta- 
ble V are not as important as the trends they reveal as 
the growth conditions are varied. 

We first consider the morphologies at high phosphine 
flow rate (R = 450). The trends in the left panel of 
Fig. 15 are similar to those seen on GaAs(OOl) (cf. Fig. 3 
of Ref. [15]), which was also grown under group- V-rich 
conditions. The direct incorporation of the adatom from 
the precursor occurs at a greater rate from above the step 
than below for both systems, but there is no asymme- 
try in the adatom step incorporation kinetics. Indeed, 
large (i.e. several tenths of electron volts) variations of 
the energy barriers for the adatom diffusion and incor- 
poration only weakly affect on the stability calculation, 



TABLE V. Frequency prefactors (u) and energy barriers (E) 
for the Arrhenius rates of the kinetic parameters in (l)-(6) 
used for the fits in Fig. 15. The parameters whose values are 
enclosed in parentheses do not affect the results in Fig. 15(a), 
even for large variations. 





R = 


150 


R = 


450 


v (s- 1 ) 


E (eV) 


v (s- 1 ) 


E (eV) 


D M 


10 14 


0.1 


10 14 


0.2 


D a 


10 12 


0.44 


10 12 


(0.6) 


K 


10 10 


0.8 


10 10 


0.8 


-1 

T 


10 10 


1.2 


10 10 


1.2 


P~M 


10 12 


0.8 


10 12 


0.8 


Pm 


10 12 


1.0 


10 12 


1.0 


Pa 


10 12 


0.75 


10 12 


(1.0) 


PI 


10 12 


0.5 


10 12 


(1.0) 



provided that the incorporation kinetics at the step edge 
are symmetric. We have enclosed the barriers for these 
processes in parentheses in Table V to indicate this fact. 
For these growth conditions, the instability of step flow 
to step bunching in InP(OOl) and GaAs(OOl) is driven 
by an effective negative step-edge barrier induced by the 
asymmetric precursor incorporation kinetics. 42 

An altogether different scenario emerges from the sta- 
bility analysis at low phosphine flow rates, as can be seen 
immediately in the right panel of Fig. 15. At lower tem- 
peratures (Tq = 520°C and 560°C) the morphology of 
all substrates shows either discrete or diffuse islands (Ta- 
ble I). Discrete islands indicate a growth regime far re- 
moved from step flow. The adatom density on the terrace 
is large, leading to a high probability of island nucleation 
and growth on the terraces. As the temperature is in- 
creased, the adatom density on the terraces decreases, 
so the growth nucleation and growth of islands decreases 
accordingly, with the islands showing some coalescence 
with steps, thereby indicating a transition regime be- 
tween discrete islanding and step flow. The separation 
of stable and unstable regions of step flow is expected to 
be qualitatively different from that obtained with higher 
phosphine flow rates. Table 5 shows that, while the asym- 
metry in the molecular step-edge incorporation rates has 
been maintained, an asymmetry in the atomic step-edge 
incorporation rates has been introduced, but in opposi- 
tion to the molecular incorporation rates. This is the 
reason for the qualitative difference in shape of the curve 
that separates the stable and unstable regions (Rk = 0). 

Our calculations show that the profile of the curve 
Rk = separating the stable and unstable regions is very 
sensitive to the atomic parameters. With decreasing tem- 
perature, the net incorporation current at steps becomes 
more negative, while that due to precursor incorporation 
becomes more positive. For certain growth conditions, 
the atomic current dominates, which suggests the onset 
of another instability through the appearance of mounds, 
a fact that has also been noted by Vladimirova et at 43 
In fact, mounds have been observed on InP(OOl) dur- 
ing mctalorganic molecular-beam 44 and chemical-beam 45 
epitaxy. The interesting point from our perspective is 
that this morphology is triggered by the low phosphine 
flow rate and can therefore be avoided. 



V. CONCLUSIONS 

Homocpitaxy on perfectly oriented substrates resulted 
in the most homogenous surfaces - stable step flow was 
observed in the majority of cases over some temperature 
range, with step-bunching composed of a maximum of 
three steps observed near 650° C. The modification of the 
V/III ratio had the most striking effect on the surface 
morphology. For R = 150, a growth temperature below 
600°C was found to be insufficient to create useful surface 
organization for any misorientation of the substrate (the 
exception being the substrate with A-steps with a mis- 
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FIG. 15. (Color online) The morphologies in Table I for different misorientation angles (toward [111]A) for R = 450 (left panel) 
and R — 150 (right panel). The contours of the decay rate in (12), which are indicted in blue, have been calculated with 
the Arrhenius parametrizations in Table V. The emboldened curves correspond to Rk = 0, which separates regions of stable 
(Rk < and (Rk > 0) step flow, and the gray regions indicate the accuracy of ±0.02° on the misorientation angles. 



orientation of 0.6°, where regular step organization was 
observed at 585°C). Temperatures above 685°C led to 
defects related to three-dimensional growth. The mid- 
range growth temperatures provided more of a distinc- 
tion between various substrates. Step-bunching was ob- 
served in most of the cases but, depending on the sur- 
face misorientation and exact growth temperature, dif- 
ferent features were formed on the surface. When exam- 
ined in detail, the overlap of crystallographic steps could 
take diverse forms, like clustering of just two monolay- 
ers (on perfectly oriented substrates) , well organized, pe- 
riodic micron-scale bunching (cliffs observed on A-type 
substrates at 655°C), or a braid-like bunching of up to 
40 monolayers (on all vicinal surfaces grown at 630° C). 
This limited the optimal growth temperature to a quite 
narrow window of about 100°C. By changing the phos- 
phine flow to R = 450, most of the samples were elim- 
inated by that restriction - stable, organized step-flow 
was observed even at 520°C and good surface morphol- 
ogy was obtained for growths conducted at temperatures 
of 720°C. 

There are several aspects of the surface morphology 
that can be addressed with modelling. The basis of our 
theoretical work 15 has been a continuum model of the 
Burton, Cabrera, and Frank 7 type. This allows an ana- 
lytic solution to be obtained and a linear stability anal- 
ysis to be carried out, but cannot address the variety 
of morphologies seen for fully developed step bunching, 
for which the solution of the two-dimensional problem 



is required. There are two approaches. One is to use a 
continuum formulation, for example, one based on the 
phase-field method. 46 This has the advantage of build- 
ing the extended terrace length scales on the substrate 
from the outset, but suffers from the difficulty of includ- 
ing the details of step kinetics while retaining a direct 
connection to the atomistic parameters. Kinetic Monte 
Carlo (KMC) simulations provide an obvious alternative 
in that the atomistic processes can be included to what- 
ever level of detail required, but the large terrace lengths 
present a computational challenge for traditional formu- 
lations. One alternative is to use a hybrid scheme that 
combines the flexibility of the KMC approach and the ex- 
tended length and time scales of continuum methods. 47 
These methods are currently under investigation. 
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